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Equilibrium properties of long-range interacting systems on lattices are investigated. There was 
a conjecture by Cannas et al. [Phys. Rev. B 61, 11521 (2000)] that the mean-field theory is exact 
for spin systems with nonadditive long-range interactions. This is called "exactness of the mean- 
field theory" . We show that the exactness of the mean-field theory holds for systems on a lattice 
with non-additive two body long-range interactions in the canonical ensemble with unfixed order 
parameters. We also show that in a canonical ensemble with fixed order parameters, exactness of the 
mean-field theory does not hold in one parameter region, which we call the "non-mean-field region." 
In the non-mean-field region, an inhomogeneous configuration appears, in contrast to the uniform 
configuration in the region where the mean-field theory holds. This inhomogeneous configuration 
is not the one given by the standard phase separation. Therefore, the mean-field picture is not 
adequate to describe these states. We discuss phase transitions between the mean-field region and 
the non-mean-field region. Exactness of the mean-field theory in spin glasses is also discussed. 



I. INTRODUCTION 

Long-range interactions cause several peculiar features: 
negative specific heat in the microcanonical ensemble [1] , 
long-lived metastable states [2], and ensemble inequiva- 
lence [3]. The nonadditivity of the interaction potential 
causes such anomalous properties. When the integral of 
an interaction potential diverges at the large distances 
(i.e., a pair interaction behaves like l/r a with a < d 
where d is the dimension of space) or when the interac- 
tion range is of the order of the system size, the system 
cannot be divided into thermodynamically independent 
subsystems. In this case, the system is said to be nonad- 
ditive. 

Until now, the statistical physics of long-range inter- 
acting systems have attracted attention [4, 5] . To under- 
stand the statistical mechanics of long-range interacting 
systems, mean-field (MF) models are employed in many 
works (see [6] and references therein). In MF models, all 
the constituents interact equally with each other regard- 
less of the distance. It is expected that the MF models 
capture some qualitative features of general long-range 
interacting systems. More strongly, some evidence of the 
exactness of the MF theory has been reported for sev- 
eral models with the power-law interaction l/r a with 
< a < d [7-12]. Here, by the phrase "exactness of 
the MF theory" we mean that the equilibrium properties 
of the system are equivalent to those of the corresponding 
MF model. In other words, the free energy of the system 
is identical to the MF free energy. Cannas et al. con- 
jectured the exactness of the MF theory for the classical 
spin systems with long-range interactions [7] and several 
studies on specific models have followed it. 

A previous work [13] demonstrated that the exactness 
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of the MF theory can be violated in a parameter region 
called the "non-MF region" for conserved systems. Ho- 
mogeneous configurations become unstable and a kind of 
phase transition should occur at the boundary between 
the MF region and the non-MF region. In the van der 
Waals limit (see below), it is well known that the in- 
stability of the homogeneous states leads to phase sep- 
aration and the configuration becomes inhomogeneous. 
These inhomogeneous states in the van der Waals limit 
can be considered as the coexistence of the two indepen- 
dent homogeneous phases, which are described by the 
MF theory. On the other hand, when the interaction is 
nonadditive, inhomogeneous states in the non-MF region 
cannot be described by the phase separation of two inde- 
pendent homogeneous phases, and they are not described 
by the MF theory. The aim of the present paper is to in- 
vestigate the nature of this phase transition between the 
MF phase and the non-MF phase. 

Although the previous work is concerned with pure 
ferromagnetic systems, we also give an extension of the 
result to the spin glass systems. 

The present work is organized as follows. Details of 
the model and setting are given in Sec. II. In Sec. Ill, 
"the exactness of the MF theory" is examined for non- 
conserved systems (systems where the magnetization is 
not conserved) and for conserved systems (systems where 
the magnetization is conserved). The exactness of the 
MF theory is always true for nonconserved systems but 
not necessarily correct for conserved systems. To demon- 
strate the above result, the Ising model with a long- 
range interaction is considered as an example in Sec. IV. 
In Sec. V, the nature of the phase transition between 
the MF phase and the non-MF phase in conserved sys- 
tems is investigated. An application of our result to spin 
glasses is discussed in Sec. VI. The summary of this pa- 
per is presented and some future problems are discussed 
in Sec. VII. 
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II. SETTING 

We consider the following Hamiltonian on a d- 
dimcnsional lattice, 



H 



T N 



N 



(1) 



Here we assumed the two-body long-range interaction. 
We impose periodic boundary conditions and interpret 
the distance between the lattice points i and j, \n — 
rj\, as the shortest distance of these lattice points in 
periodic boundary conditions. The lattice interval is set 
to unity. The parameter J > is a coupling strength 
and K(Ti — rj) is the interaction potential between the 
sites i and j. When ^2 r .^. K(ri) diverges, the system 
is nonadditive and there exists no thermodynamic limit 
in the usual sense. To avoid this difficulty, we normalize 
the interaction potential as 



(2) 



This is called Kac's prescription in the literature. The 
"spin" variable at is arbitrary as long as it is finite; in 
the Ising model cr; = ±1, in the classical XY model 
<7i = (cos 9, sin 9), and in the g-state Potts model Ui = e a , 
a = 1, 2, . . . , q, where e a ■ = S a b, and so on. In this 
paper, we treat the one-component spin variable to make 
the presentation simple, but the generalization to multi- 
component spin variables is straightforward. 

As the simplest long-range interacting model, the 
infinite-range model exists, 



J 



N 



N 



^MF = - #E CTi ' 



(3) 



for which it is known that the MF theory is exactly ap- 
plicable. Hereafter we call this model the "MF model." 

In present paper, we consider the following two types 
of long-range interactions: the power-law interaction 



K(r) oc — , < a < d, 

and the Kac potential [14] 

K(r) cx ~/ d 4>(jr). 



(4) 



(5) 



Here, <fi(x) is assumed to be non-negative <fi(x) > and 
integrable / d d x(j)(x) < +oo. Moreover, we assume that 
there is a positive and decreasing function ip(x) such that 



/ 



\V<P(x)\ <V(*), 
d d xijj(x) < +oo. 



(6) 



This assumption is necessary to justify the coarse- 
graining of the Hamiltonian discussed later. A typical 



example of the Kac potential is the exponential form, 
K{vi — r.j) oc 'y d e~' r \ ri ~ r i L In this case, 4>{x) = ip{x) — 
e~ x . 

We will take the limit 7 — >• in the Kac potential. In 
this paper, two limiting procedures are considered: the 
van der Waals limit 7^0 after L — > 00 [15] and the 
long-range limit 7^0 with jL = const. The former limit 
corresponds to the situation where the interaction range 
7 _1 is much longer than the microscopic length scale (the 
lattice interval) but much shorter than the system size L. 
In this case, the system is additive and it does not show 
anomalous behavior like the ensemble inequivalence. The 
latter limit corresponds to the situation where the inter- 
action range is comparable to the system size. In this 
case, the system has no additivity. These two limits give 
different behavior in general. 



III. EXACTNESS OF THE MF THEORY 

A condition of the exactness of the MF theory has been 
reported briefly [13]. In this section, we give the detailed 
explanation of this property. 

In the long-range interacting systems, it is expected 
that only long wavelength modes play important roles for 
macroscopic behavior. In fact, it is possible to perform 
coarse graining exactly for long-range interacting models. 

Now we explain what coarse graining is. Let us divide 
the lattice system into blocks of the linear dimension I. 
The number of blocks is {L/l) d and each block has l d 
sites. We introduce the local coarse-grained variable m p 
as 



m, 



(7) 



ieB p 



in each block B 
the limit L 



p , where p = 1,2,..., (L/l) d . We take 
00, I — > 00 with l/L — > (continuous 
limit). This strategy is the same as the procedure in 
the paper by Barre et al. [10]. We define the position 
x p = r p/L, where r p is the central position of a block 
Bp \p = 1, . . . , (L/l) d \. We also define m(x p ) = m p . For 
long-range interacting models, as shown in Appendix A, 
the Hamiltonian is expressed only by m(x) in the ther- 
modynamic limit: 



H = U[m(x)] + o(N), 



(8) 



where 



H[m{x)] = - [ d d x [ d d yU(x - y)m(x)m(y) 
2 Jc d J c d 



(9) 



NH / d a xm(x). 
Jc d 



Here, the scaled interaction potential U (x) is given by 
U(x) = lim L d K(Lx). (10) 
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The integrations in Eq. (9) are performed over a d- 
dimensional unit cube Cd, namely x,y G [0, l] d . Kac's 
prescription (2) implies 



U(x)d d x = 1. 



(11) 



In the power-law potential and in the Kac potential 
with the long-range limit, U(x) = K\jx a and U{x) = 
K27q0(7o:e), respectively. Here, ki and k-i are normal- 
ization constants determined by Eq. (11) and 70 = jL is 
a constant in the long-range limit. In the Kac potential 
with the van der Waals limit, U(x) = 8{x). 

Performing the Fourier expansion in Eq. (9), we obtain 
the following expression: 



H = -—^y^Un\m n \ 2 



NHm a , 



where 



m(x) — m n e 



2Ttin-x 



u n = 



/ d d xU{x) cos(27rn • x) 
Jc d 



(12) 

(13) 
(14) 



We call {U n } interaction eigenvalues. Interaction eigen- 
values of n ^ are less than or equal to unity U n < 1, 
because 

|Z7 n | = / d d xU(x) cos(2?rn • x) 
Jc d 

< / d d x\U(x) \ ■ |cos(27rn-a;)| 
Jc d 

< I d d xU(x) = 1. (15) 
Jc d 

From now on, we consider the generalized free energy 
F(m,T, H) which is defined as 



exp[-/3F(m,T, H)} 



E 



(16) 



[{<Ti}\m =m] 



where the summation is taken over the configurations 
with a fixed value of the total magnetization mo = to. 
The temperature is T — 1/(3. We can separate the long- 
wavelength modes from the short ones by coarse graining: 



E 

[{ff,}|m =m] 



Vm{x) ^2 ■ (I 7 ) 

1,1 "' i&i} with fixed m(x) 



The summation with the fixed m(x) is expressed as 

1 = exp ( / S{m(x))d d x j , (18) 

{o-i} with fixed m(x) \ J C d J 

where S(m) is the entropy, 

S(m) = ln(the number of states 

with the fixed magnetization to). (19) 




CE{g(x)} 



FIG. 1: An illustrative example of the convex envelope of the 
function g(x). The solid line denotes g(x) and the dashed line 
denotes CE{g(x)}. 



From Eqs. (8), (17), and (18), we obtain 
exp[-l3F(m,T,H)}= Vm{x) 

J rho— m 

x exp -P (ii\m{x)} -T J S(m(x))d d x^ . (20) 
By using the saddle-point method, we have 



F(m,T,H)— min 

{m(x)\fho—m} 



H[m(x)] - T I S{m{x))d d x 
c d 

(21) 



min T({m(x)},T,H) 

\m{x)\mQ—7n\ 



The functional T is called the free energy functional here- 
after. See Appendix B for the rigorous justification of 
using the saddle-point method. 

From Eq. (21) we obtain the upper bound 



F(m, T, H) < F[{m{x) = to}, T, H] 
= %mf — TS(m) 
= FuF(m,T,H). 

The lower bound 

F(m,T,H) > F MF (m,T, H)-U max AF MF (m,T/U n 



(22) 



(23) 



is obtained by replacing all the U n with n ^ by U„ 
which is defined as 



max Un 



The function AFmf is defined as 

Ai^MF = -Fmf — CE{Fmf}, 



(24) 



(25) 



where CE means the convex envelope. The convex enve- 
lope of a function g(x) is defined as the maximum convex 
function not exceeding g{x) (see Fig. 1). The derivation 
of the lower bound (23) is given in Appendix C. 

Collecting the upper and the lower bounds, we have 

Fmf(to, T, H) - U ma , x AF M F(m, T cSl H) 
<F(m,T,H)< F MF (m,T,H), (26) 
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where T oS = T/U max . 

An inequality necessary to prove Eq. (23) is Eq. (C3) 
which corresponds to the replacement of U n — > f/ max 
for all n ^ 0. In the Kac potential with the van der 
Waals limit, U n = 1 for all n. Therefore, the equality 
F(m, T, H) = CE{F MF (m,T,H)} holds for the Kac po- 
tential with the van der Waals limit. It is a well known 
fact that the MF theory with the Maxwell's equal area 
law is justified in the van der Waals limit [15]. The 
Maxwell's equal area law is equivalent to the replace- 
ment of the MF free energy by its convex envelope. The 
replacement of the free energy by its convex envelope in- 
dicates the occurrence of the phase separation. 

Next, we consider the local stability of the uni- 
form configuration m(x) = m for all x. If the uni- 
form configuration gives the local maximum of the 
free energy functional, there are configurations which 
have lower free energy than the MF. Therefore, in 
this case, F(m,T, H) < Fmf{iti,T, H) holds instead of 
F(m, T, H) < F M F(m,T,H). From this local stability 
analysis, we have 

F(m, T, H) < F MF (m,T,H) 

d 2 

for ^F MF (m,T cS ,H) < 0. (27) 

We have derived necessary inequalities for the gener- 
alized free energy with this. According to Eqs. (26) and 
(27), the parameter region (m,T,H) is classified to the 
following three regions: 

Region A : the region where AFmf (w,T e ff,_ff) = 0. In 
this region, the MF model gives the exact free en- 
ergy, F(m, T, H) = F MF (m, T, H). 

Region B : the region where AFmf (wi, T e s, H) > and 
~M^?Fmf {m,T c g, H) > 0. In this region, it is not 
sure whether the MF model is exact, F(m, T, H) < 
Fyi F {m,T,H). However, the homogeneous states 
determined by the MF theory are locally stable. 

Region C : the region where g^s-FiyrF (w, T e gp, H) < 0. 
In this region, the MF model cannot describe 
the long-range interacting systems, F(m, T, H) ^ 
Fw£p(m,T,H). In this region, the homogeneous 
states are not even locally stable. 

Notice that this classification is determined only by 
Fmf and J7 max . Hence, we can specify these three regions 
concretely for individual models by analyzing only the 
MF models. 

In region A and a part of B, F(m,T,H) = 
F^i F {m,T,H) holds. This region is called the MF re- 
gion. On the other hand, in region C and the other part 
of B, F(m, T, H) ^ F MF {m,T,H). This region is called 
the non-MF region. 

In region B, we cannot say whether a point (m, T, H) 
belongs to a MF or non-MF region. It depends on the 
type of the interaction, the value of the temperature, and 



the specific model (what <7j is). However, the homoge- 
neous state described by the MF theory is locally stable. 

In the non-MF regin, some inhomogeneous modes must 
develop. Therefore, by observing the equilibrium config- 
uration of the system with the conserved magnetization 
m, we can know whether this point belongs to the MF 
or the non-MF region. Namely, if the cluster appears in 
equilibrium, this point (m, T, H) belongs to the non-MF 
region; on the other hand, if the system is uniform, this 
point belongs to the MF region. 

Here we discuss the exactness of the MF theory 
based on the derived inequalities. In nonconserved 
systems, the equilibrium magnetization rn C q is deter- 
mined by min m [F(m,T, H)] — F(m oq ,T, H). Because 
the equilibrium state belongs to the MF region where 
AJ 7 MF(w eq , T e g) = 0, it is concluded that exactness of 
the MF theory is true at any temperature in nonconserved 
systems. 

In contrast, in conserved systems, the generalized free 
energy itself is the equilibrium free energy and the de- 
rived inequalities mean that the equilibrium property of 
the long-range interacting system is exactly the same as 
that of the corresponding MF model in the MF region. 
On the other hand, they are not the same in the non-MF 
region. As discussed above, the inhomogeneity appears 
in the non-MF region. The clustering phenomena can- 
not be described by the MF model with the help of the 
standard phase-separation argument. Therefore, we con- 
clude that exactness of the MF theory is violated in con- 
served systems. We investigate the clustering phenomena 
in Sec. V. 



IV. EXAMPLE: LONG-RANGE INTERACTING 
ISING MODEL 



T>J T<J 




FIG. 2: The MF free energy of the Ising model (the solid line) 
and its convex envelope (the dashed line). The left figure is 
for the case T > J. As the MF free energy is convex for all 
m, it is equivalent to its convex envelope. The right figure is 
for the case T < J. In this case, the MF free energy is not 
convex, and the flat region appears in the convex envelope. 



To make the statement clear, let us pick the long-range 
Ising model as a simple example. We consider the Hamil- 
tonian (1) with <x; = ±1. The corresponding MF model 



■5 



is 

(28) 

We can calculate the MF free energy, which is given by 

/mf = - ^m 2 - Hm 

„ / 1 + m , 1 + m 1 — m , 1 — m\ . , 

+t (— in — +— in — )■ < 29 > 

In Fig. 2, the MF free energy and its convex envelope 
at H = are depicted. For T < T c = J, it becomes a 
nonconvex function of to. 

When the "effective temperature" T e g = T/U max is 
above the critical temperature, T c = J, the relation 
F MF {m,T eS ,H) = CE{F MF (m,T eS ,H)} holds for any 
to. In this case, Eq. (26) leads to 

F(m,T,H) = F MF (to, T, H) (30) 

for any to. On the other hand, when T e g < T c , there is a 
region where the MF free energy is not convex, and from 
Eq. (26) we can conclude that there is a region of to such 
that F(m, T, H) ^ F MF (to, T, H). 

A schematic picture of the boundaries of three regions 
is depicted in Fig. 3. The upper line is the MF free energy 
at the temperature T e g = T/£/ max and the lower line is 
that at the genuine temperature T. Figure 3(a) describes 
the case of T e g > T c > T. In this case, Fmf(to, T e s, H) is 
a convex function of to. Therefore, A^mf(w, T c g, H ) = 
for any to and H. Figure 3(b) describes the case of 
T c > T c ff > T. In this case, Fmf(wi, T e s, H) is not con- 
vex and deviation from the convex envelope (dashed line) 
appears. 

In the Ising model with a long-range interaction, we 
can give the expression of regions A, B, and C explicitly. 
In region A, AFmf("i, Toff, H ) = holds. Because AFmf 
is independent of H, the region A is given by 

\m\>\m cq (T cS ,H = 0)\, (31) 

where m cq (T,H) is the magnetization in equilibrium 
which is given by the self-consistent equation, 

TO cq (T, H) = tanh [/3(Jm eq (T, H) + H)] . 

Similarly, region B is given by 

|m sp (T eff )| < |to| < |TO cq (T cff ,0)|, (32) 

where m sp (T) is the spinodal point that is the metasta- 
bility limit in the nonconserved systems, 

m sp (T) = ± \j l -j: J - ( 33 ) 
Finally, region C is given by 

\m\ < K P (T cff )|. (34) 



After all, m AB = m cq (T off ,0) and m B c = m sp (T cS ) in 
Fig. 3. 

In conserved systems, a typical spin configuration is 
homogeneous at least in region A [see Fig. 4(a)] but is 
inhomogeneous in region C [see Fig. 4(b)], as predicted 
in the previous section. 

In nonconserved systems, we confirm the MF model 
is exact. We depict spontaneous magnetizations for the 
two-dimensional Ising model with a long-range interac- 
tion (a = 1) in Fig. 5. As predicted in the previous 
section, it agrees with the MF result determined by solv- 
ing the self-consistent equation, to = tanh/^Jm. More- 
over, it turns out that not only the equilibrium states 
but also the metastable states of the MF model given 
by the local minimum of the free energy are maintained, 
because the local minimum is located in the range of 
|w| > |to sp (T)|(> TO sp (r c ff)|); namely, it is located in re- 
gion A or B. In region A or B, homogeneous states are lo- 
cally stable against the inhomogeneous fluctuations, and 
therefore the metastability defined in the MF model is 
not lost. If they belonged to region C, they would have 
instability against inhomogeneous fluctuations and lose 
their local stability. 



V. PHASE TRANSITION BETWEEN MF 
PHASE AND NON-MF PHASE 

A kind of phase transition takes place between MF 
and non-MF regions in conserved systems. We demon- 
strate a typical spin configuration of the Ising model with 
1/r type long-range interaction in Fig. 4. In Fig. 4 (a), 
we depict a configuration at T/J = 0.28 and to = 0.6 
which is in the MF region and we find a homogeneous 
configuration. In Fig. 4 (b), we depict a configuration at 
T/J = 0.18 and to = 0.6 which is in the non-MF region, 
where we find a clustered configuration. We confirmed 
that there are two parameter regions, that is, the MF 
region and the non-MF region. 

The next problem is to understand mechanism of the 
phase transition between the MF phase and the non-MF 
phase in conserved systems. The Fourier modes with 
large interaction eigenvalues will play a significant role 
in the phase transition. 

First consider the Kac potential with the van der Waals 
limit. In this case, all the interaction eigenvalues are 
unity, U n = 1; all the Fourier modes contribute to the 
phase transition. This fact indicates that the standard 
phase separation occurs and the system can be divided 
into two subsystems with different phases. 

On the other hand, in the power-law interactions or 
the Kac potential with the long-range limit, the spec- 
trum of the interaction eigenvalues is discrete; only the 
small number of Fourier modes which have the maximum 
interaction eigenvalue will contribute to the phase tran- 
sition. Therefore, it is expected that the phase transition 
is understood by the Landau expansion of the free energy 
functional by such Fourier modes. It is also expected that 
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m~AB 'mac 

FIG. 3: (Color online) An illustrative explanation of the relation of regions A, B, and C together with the MF free energy 
in the Ising model. The upper line is the MF free energy at the temperature T e a and the lower line is that at the genuine 
temperature T. (a) The case of T cS >T C >T. (b) The case of T c > T cff > T. 




FIG. 4: Typical equilibrium snapshots in the two-dimensional (2D) lattice gas model (conserved Ising model) with 1/r-type 
long-range interaction. Periodic boundary conditions are imposed. Black points represent occupied sites, (a) Region A. The 
parameters are set to be m — 0.6, T = 0.28, J = 1, and L — 100. (b) Region C. The parameters are set to be m = 0.6, 
T = 0.18, J — 1, and L = 100. 



MF solution 




0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 

T 

FIG. 5: Magnetizations against temperatures. The points 
are numerical results of the 2D long-range Ising model with 
K(r) oc 1/r and L — 128 calculated by the Monte Carlo 
method. Error bars are smaller than the symbol size. The 
solid line is the MF solution. They agree very well. 



this phase transition is qualitatively different from that 
in the Kac potential with the van der Waals limit, which 
is well described by the notion of the phase separation. 



We focus on the modes with the maximum interaction 
eigenvalue U max , and the set of these modes are denoted 
by-M, 

M = {n:n€Z d ,U n = U max }. 

We regard rh n with tieMas order parameters, and the 
other modes are determined by the condition 

8H{m{x)},T,H) =Q 

which is the equilibrium condition under the fixed order 
parameters. 

In many cases, including the power-law interactions, 
U n tends to decrease as the length of the wave number 
vector | to | increases. Therefore, we assume that 

M = {±e a }, a=l,2,...,d (36) 

where e Q is a unit vector and e a -e.b = 5 a b- In this section, 
we assume that the interaction is isotropic. Hence, the 
order parameters m± ea do not depend on the direction 
a, so we put <j> = m± ea and we regard 4> as an order 
parameter. Let us expand the free energy functional by 
{m± ea } with the help of Eq. (35). It is noticed that 

m n = O(|0|l«il+KI+-"+M), (37) 
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which is derived from Eqs. (35) and (36). Up to the some calculations: 
fourth order of \<p\, we obtain the following expansion by 



( [d (T^) 2 

Fmf + N j-^JC/i + Ts^M 2 + [- A j U2 + t s{ 2) + d ( d 



1) 



d(d- 



+o(H 6 ). 

(38) 



r 



We defined = d k s(m)/dm k and 



Ui = U ±ea = U n 



U n 
U 2 



Ue a ±e b 
U±2e a - 



(39) 



Here we rewrite Eq. (38) as 



N 



/MF+a|^| 2 + 6|0| 4 + O(|0| 6 ). 



(40) 



It is noted that a = corresponds to 
<9 2 -Fmf(w, T e ff, H)/dm 2 = 0, that is, the boundary 
between regions B and C. It is evident that the first 
order phase transition occurs if b < when a = 0. It 
is reasonably expected that the second-order transition 
occurs at the boundary of regions B and C if b > 
when a = 0, though we do not show it strictly. If we put 
a = 0, then 



dJU x 



s (4) 

(M-l)^ 



4(d-l) 



+ 



,(3) 



=6o- 



(41) 



From the above discussion, the transition is of the first 
order when b < and of the second order when bo > 0. 

Let us demonstrate the above result in the Ising model 
(<7j = ±1). In this case, the entropy is given by 



s(m) = — 
In this model 



1 + m , 1 + m 
In 



m , 1 
— In — 



TO 



(42) 



6 oc(2d- 1)(1 + 3m 2 ) 



2to^ 



■4(d- 1)- 



=(2d- 1)(1 + 3m 2 ) - 2to 2 /sT. 
If we put 6o = 0, 

2d - 1 



to 2 = 



2A"-3(2d-l) 



= to 2 



(43) 



(44) 



The first-order transition occurs when \m\ > \m c \ and the 
second-order transition occurs when \m\ < \m c \. When 
Un, U 2 — > Ui, then b is negative for any m ^ 0. There- 
fore, the transition is first order except for to = in 
the case where the spectrum of interaction eigenvalues 
is almost continuous. On the other hand, if we put 
U\\ = U 2 = 0, bo is positive and the transition is second 
order for any to. Thus, we have found that the long- 
range nature of the interaction (the discreteness of the 
spectrum of interaction eigenvalues) enhances the second 
order phase transition between homogeneous and inho- 
mogeneous phases. 

As an example, let us consider the Ising model on 
the two-dimensional square lattice with the interaction 
K{r) - 1/r. In this case, U x w 0.310, U n w 0.207 
and U 2 « 0.132. These parameters imply m c m 0.402. 
We demonstrate the result in Fig. 6. In this figure, we 
plot the average values of the order parameter \<f>\ in the 
Monte Carlo dynamics under a temperature sweep over 
the range T = 0.1 J to 0.4 J and T = 0.4 J to 0.1J. The 
transition is continuous when m = 0.3 but it is discon- 
tinuous when to = 0.5 and to = 0.6. Indeed, hysteresis 
loops appear for to = 0.5 and 0.6 as in Fig. 6. 

Thus, the phase transition between the MF and the 
non-MF phases does not necessarily occur just at a = 0, 
the boundary of the regions B and C. There are situations 
that the homogeneous states are locally stable but not 
globally stable. This shows the relevance of considering 
the region B. 



VI. APPLICATION TO SPIN GLASSES 

In this section, we discuss an application of the result 
discussed in Sec. Ill to spin glass systems, whose Hamil- 
tonian is 



T~L — — ^""^ Kjj JijUjO j , 



(45) 



i<3 



where Kij denotes the two-body interaction potential 
and Jij is assumed to be a Gaussian random variable 
whose probability distribution P(Jij) is given by 

1 ( 1 



P(Jii) 



(46) 



8 




FIG. 6: Monte Carlo results of average values of the order 
parameter \<p\ with the temperature sweep from T = 0.1J to 
0.4J (1-h) and from 0.4J to 0.1J (h-1). The open circle is the 
data for m — 0.3 and 1-h, and the closed circle is for m = 0.3 
and h-1. The open square is for m = 0.5 and 1-h, and the 
closed square is for m = 0.5 and h-1. The open triangle is 
for m = 0.6 and 1-h, and the closed triangle is for m = 0.6 
and h-1. The system size is L = 60 for m = 0.3 and L = 100 
otherwise. 



The long-range spin glass systems have been extensively 
studied recently in order to extend our knowledge of spin 
glasses in finite dimensions [16-19]. In these works, the 
universality class of the spin glass systems with 1 / r a -type 



interaction has been studied. In one dimension, it has 
been revealed that the transition is in the MF universality 
class in the case of 1/2 < a < 2/3 and the non-MF 
universality class in the case of 2/3 < a < 1. In this 
paper, we focus on the nonadditive regime, a < d/2, and 
show that the system is fully identical to the MF model 
(Shcrrington-Kirkpatrick model [20]) at any temperature, 
not only the critical exponents. 

The free energy of this system is expressed as 



lim 

N-toc 



1 

7V/3 



(InZ). 



(47) 



Here the angular bracket denotes the average over the 
random interactions { Jy } and Z is the partition function, 
Z = Tre~^ w . In this section, the symbol Tr refers to the 
summation over all the microscopic configurations {<Tj}. 

In order to examine the free energy, we apply the 
replica method, which is a nonrigorous but successful 
technique [21]. We can express the free energy as fol- 
lows: 



1 ( 7 n \ 
f(B) = - lim lim J '- 



1 



(48) 



In the replica method, first we calculate (Z n ) for integer 
n, then extrapolate it to noninteger n and take the limit 
n — > 0. 

Let us calculate (Z n ) for integer n. The quantity (Z n ) 
is expressed as 



(Z n ) = Tr J m 



dJi 



exp 



— V J 2 

2J 2 J v 



r 



a—l i<j 



(49) 



By integrating out over {J^}, we obtain 

(Z n ) =Trexp ( i/3 2 J 2 £ E ) 

y a, b i<j J 

If we define the vector S by 

(Si) n ( a -l)+b = ^i^i, 

then we obtain the formal expression, 
If we define the following quantities, 



we obtain 



(50) 



(51) 



(52) 



(53) 



(Z n ) = Trexp 



(54) 



In this form, we can apply the argument studied in the 
present paper on the exactness of the MF theory for 
pure ferromagnetic systems to spin glasses. Namely, it 
is straightforward to show that MF theory is exact as 
long as the interaction Kij = Kfj is long range. If we 
assume the power- law interaction Kij = C/rfj, then the 
MF model is exact when a < d/2 instead of a < d. The 
constant C is determined by the normalization condition 



r 2 

V— = 1 



(55) 



which is the correspondence of Eq. (2). (Note that in the 
infinite-range case, C = 1/ y/N.) 
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From the above argument, we can conclude that (Z n ) 
in the system (45) is exactly equal to that in the MF 
model at least for any integer n. Therefore, the two free 
energies calculated by the replica method may be exactly 
equal. Unfortunately, it is not certain if this result can be 
proved without using the replica method. However, from 
the above result, it is reasonably expected that the true 
free energy of the system (45) is also exactly equal to that 
of the corresponding infinite-range model. If the equilib- 
rium properties of the system (45) with = C/rfj and 
a < d/2 in contact with a thermal reservoir at a tempera- 
ture T were different from those of the corresponding MF 
spin glasses, then it would imply that the replica method 
is not exact at least for the long-range interacting spin 
glasses. However, we have no reason to expect that the 
replica method does not work in long-range spin glasses, 
and we conclude that spin glass systems with the Hamil- 
tonian (45) and a < d/2 exhibits behavior identical with 
the corresponding MF models. 



VII. SUMMARY AND DISCUSSION 

We investigated the exactness of the MF theory in sys- 
tems with nonadditive long-range interactions (power- 
law potential or the Kac potential with the long-range 
limit) and additive long-range interactions (the Kac po- 
tential with the van der Waals limit) in a unified way. We 
showed that the exactness of the MF theory is always 
valid for the nonconserved systems, while for the con- 
served systems there exists the parameter region where 
the exactness of the MF theory is violated. As an ap- 
plication of our result, we considered the spin glass sys- 
tem and revealed that the exactness of the MF theory 
is valid also for long-range interacting spin glass systems 
within the treatment by the replica method, as long as 
the square of the interaction potential is nonadditive. 

We examined the nature of the phase transition be- 
tween the MF phase and the non-MF phase by the Lan- 
dau expansion of the free energy functional. We pointed 
out that except for the van der Waals limit, inhomoge- 
neous states observed in the non-MF region are quite dif- 
ferent from those created by the phase separation. This 
aspect is reflected to the fact that only a small number of 
Fourier modes are important in the inhomogeneity. It is 
indicated by the discrete spectrum of interaction eigen- 
values U n . On the other hand, in the van der Waals limit, 
all the interaction eigenvalues are degenerate, U n = 1, 
and the standard phase separation occurs. It will be in- 
teresting to study dynamical nature associated with the 
phase transition between the MF phase and the non-MF 
phase, because it is known that long-range interacting 
systems exhibit peculiar features also in dynamics. 

The difference between conserved and nonconserved 
systems is a consequence of the violation of ensemble 
equivalence. Whether the exactness of the MF theory 
holds or not depends on the specific ensemble. A natural 
question is whether the MF theory is exact for a mi- 



crocanonical ensemble. It was investigated in Ref. [10], 
where it turned out that the exactness of the MF theory 
is valid for the long-range Ising model in the microcanon- 
ical ensemble. Here, it should be noted that in the MF 
Ising model, the canonical ensemble and the microcanon- 
ical ensemble are equivalent. We can show from the re- 
sult presented in Sec. Ill that the exactness of the MF 
theory for long-range interacting systems is valid also in 
the microcanonical ensemble if the microcanonical and 
canonical ensembles are equivalent in the corresponding 
MF model. When the two ensembles are not equivalent 
in the MF model, however, it is not obvious whether the 
exactness of the MF theory holds in the microcanonical 
ensemble. This issue will be investigated elsewhere. 

In non-conserved systems, the MF theory is always ex- 
act in equilibrium. However, in the out-of-equilibrium 
situations, the inhomogeneity due to the long-range in- 
teractions may appear. For example, let us consider the 
relaxation from the metastable states with the uniform 
magnetization profile. The MF metastable states as lo- 
cal minimum of the free energy are also remained in the 
general long-range interacting systems, as pointed out in 
Sec. IV. However, at low temperatures these metastable 
states belong to region B. If the metastable states belong 
to the non-MF region, these metastable states may relax 
to equilibrium by appearing the temporal inhomogeneity 
Hence, in these situations, the relaxation to equilibrium 
will be different from that observed in the MF models. 

Recently, it was reported that a model of spin-crossover 
materials has an effective long-range interaction among 
molecules [22]. In this model, although the Hamiltonian 
has only short-range interactions, effective long-range in- 
teractions among molecules appear due to the lattice dis- 
tortion by the difference of the molecular size depending 
on molecular states. In such systems, the results of the 
MF model including the MF spinodal are not artifacts 
but physically relevant results. 

In small systems, the range of the interaction can be 
of the order of the system size. Indeed, negative heat 
capacity has been observed experimentally in small sys- 
tems [23]. For example, the phenomenon of the super- 
radiance originates from the effective long-range inter- 
action among two-level atoms mediated by the coupling 
with a cavity mode [24]. The similarity between long- 
range interacting macroscopic systems and small systems 
should be discussed in the future. 

In this way, the nature of long-range interacting sys- 
tems may be widely observed, and it will become more 
important to study it from a general point of view. 
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Appendix A: Justification of the coarse-graining 

We justify the coarse graining (9) in this appendix. 
The coarse-grained Hamiltonian with a finite system size 
N = L d and a finite number of blocks il = (L/l) d is given 
by 



where 



J- 

2 2-1 U pq ! " I*'"': 



-Hl d J2 m P> ( A1 ) 



ieB p jeB q 



(A2) 



and m p is the average global variable of the block B p 
defined by Eq. (7). Notice that H(N,Q) approaches the 
coarse-grained Hamiltonian (9) when the limit N — > oo, 
£1 — > oo with fl/N — > is taken. 

Our aim is to prove that there exists a sequence f2(iV) 
such that 



N 

and satisfies 



hm ^- sup \n-n{N,n(N))\ = o, 



N 



' lim Q(N) = oo, 

N— ¥00 



Q(N) 



K AT->-oo N 



lim 



0. 



(A3) 



(A4) 



When we consider the Kac potential with the van 
der Waals limit, we replace lim7v->oo and fi(iV) by 
lim 7 _>o limAr^oo and Cl(N, 7), respectively. 

The following derivation is a strightforward generaliza- 
tion of the strategy of the paper by Barre el al. [10] , which 
they proved only for the one-dimensional Ising model 
with the power-law interaction. 

First, we express the coarse-grained Hamiltonian in 
terms of the microscopic variables {<X;}. Then 



U(N,{1) 



J 
2 



i&B p ,jeB q 



(A5) 

where I is the linear dimension of a block. We thereby 
have 



J 



< 



2Nl 2d 

J_ 

2NI 



E E E [K{ni)-K{r u )] 

v,q ieB p ,j£B q keB p ,ieBg 

3d E E M E mr^-Kiru, 



p,q iGB p ,j£B q 



k^LBp ,l£LB„ 



(A6) 



























B q 



























B p 

p 6 dq 



FIG. 7: An illustrative explanation of the set dq for d — 2. 
The central block is B q and the surrounding gray blocks are 
in dq. 



Assuming |ct;<7j| < C, we have 



1 



< 



CJ 
2Nl 2d 



EE E mr^-Kiru) 



o,q i£B p ,j£B q k£B p ,l£B q 



(A7) 



We define D pq as the shortest distance between the two 



blocks Bp and B q . Namely, 

D pq = min _ \r k - rj|. 



k^Bp ,l^B q 



(A8) 



Moreover, the set dq is defined as the set of p such 
that D pq < I (see Fig. 7). We divide the summation 
of Eq. (A7) into two terms: the term with p £ dq and 
the term with p ^ dq. Then we obtain 

jj\H-*\ 



<^EE E E m rij ) - K { r kl )\ 

q=l p€dq iGB p ,jeB q kGB p ,lGB q 

+ ^EE E E l*(r«)-*(r«)| 
q=i p<£dqieB p jeB q keB p deB q 

= A l +A 2 . (A9) 

From now on, we prove that there exists a sequence 
Q(N) such that liniAr^oo A\ — liiriAr^oo A2 — both for 
the power-law potential and the Kac potential. 

Upper bound of A\ 

If we define Kq as the maximum value of \K (ry)|, then 
\K(rn) - K{r M )\ < \K{ rij )\ + \K(r kl )\ < 2K Q . (A10) 
Therefore, we obtain 



A x < 



2AZ 2d ^ ^ 

9=1 pGa<j 



P d • 2if = 



2d n 

v E 2. 

(All) 
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where we used X^es 1 = ■ Because the number of 
blocks which belong to dq is determined by the spatial 
dimension d, we can write 



E 1 = n <i- 

p£dq 

Then we obtain the upper bound of A\: 

CJK Q l 2d n 
Ai < — Yln d . 

As n = (L/l) d and N = L d , 

A l < Cn d K l d . 



(A12) 



(A13) 



(A14) 



Remember that Kq depends on the system size L or the 
interaction length 7, because we normalize the interac- 
tion potential so that 



E = L 



(A15) 



In the case of the Kac potential, K(r) ~ 'j d (f>(jr). 
Here, the symbol ~ means equal except for a nonessential 
factor independent of I, L, and 7. Therefore, the symbol 
~ does not imply any approximations (we use the symbol 
< similarly) . Thus K ~ j d and 



K l d ~ {ll) d - 



(A16) 



In the van der Waals limit, we take 7 — > after L — > 00. 
Hence if we take the limit I — > 00 o/ter the limit 7—^0, 
then Ai — > and the conditions (A4) are satisfied. 

In the long-range limit, where L — > 00 with 7L = 
const, is taken, if we take I — > 00 after L — > 00, then 
Ai — >■ and the conditions (A4) are also satisfied. 

In the case of the power-law potential, K(r) ~ 
L a - d /r a and K ~ L a - d for < a < d. Therefore, 



A 1 



K Q l d 



I 



L l-a/d 



(A17) 



After all, A\ — > when we take the limit of Z — > 00 after 
L — > 00. When a > d, the function l(L) which satisfies 
Eq. (A4) does not exist, and the coarse graining cannot 
be performed exactly. 



Upper bound of A2 

As p £ dq, D pq ^ and the following inequality is 
satisfied for n, £ B p and rj,ri £ B q : 

\rij -r k i\ = \{n - r k ) - (rj - n)\ 

< \n - rk\ + \rj - n\ 
= r ik + rji 

< l^fdl = Dl (A18) 



From the mean-value theorem, there exists < 9 < 1 
such that 

K(r i3 ) - K{r kl ) = [VK{r l3 + 9{r kl - r tj ))] ■ {r kl - r y ). 

(A19) 

Combining the above two relations, we have 

\K( rij ) - K(r kl )\ < Dl\WK( rij + 6(r kl - r«))|. (A20) 
Here, for p ^ q, 

\VK(r tJ +6{r u -r tf ))| 

(A21) 



Using these inequalities, we can evaluate the upper bound 
ofA 2 : 



CDJl 2d+1 A ^ . 

a 2 < — EE^- 

3=1 pfdq 



2N 



(A22) 



As the blocks are aligned in the (i-dimensional space, 
each block can be labeled by the <i-dimensional vector 

n = (ni,n 2 , • • • ,n d ), 

where rij is an integer and —L/21 < rii < L/21. We fix 
the block B q at n — 0, and we determine ft so that the 
configuration r p of the center of a block B p may be given 
by r p = In. In this case, there is a constant < a < 1 
such that D pq > al\ri\ for all p with B p £ dq (see Fig. 8). 
When we assume periodic boundary conditions, [25] the 
summation can be written as follows: 

n L/21 L/21 L/21 

EE^<« E E ••• E K ko 

q=lp</dq n 1 = - L/21 n 2 = - L/21 n d =-L/2l 

(A23) 



where the index p corresponds to the vector n and we 
write K' pq as KL Q as well as we defined KL = for 
n £ dq. 

Thus, from (A22), the upper bound of A 2 is given by 



^<^ +1 E4o- 



(A24) 



Next, we consider each case of the interaction forms and 
evaluate the upper bound of A 2 . 

In the Kac potential, K{r) ~ 7 <Z </K7 r ') an d therefore 

K'(r) =-f d+1 4>'(-/r). (A25) 

From the assumption of the Kac potential (6), 

K'nfi = l d+1 Tp(Dn,o) < J d+1 i>(al\n\). (A26) 
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n = (ni,n 2 ) 



(-2,2) 


(-1,2) 


(0,2) 


(1,2) 


(2,2) 


(-2,1) 


(-1,1) 


(0,1) 


(1,1) 


(2,1) 


(-2,0) 


(-1,0) 


B q 


(1,0) 


(2,0) 


(-2,-1) 


(-1,-1) 


(0,-1) 


(1,-D 


(2,-1) 


(-2,-2) 


(-1,-2) 


(0,-2) 


(1,-2) 


(2,-2) 



(1,2) 



B r 



FIG. 8: An illustrative explanation of n for d — 2. The 
distance between the central positions of and B q is given 
by l\n\. Moreover, for all n, Bft dq, there exists a constant 
< a < 1 such that the shortest distance between B^ and 
B q = Bo is restricted by Dh,o > al\n\. For d — 2, this 
constant is a = 1/yfb (see the right figure). 



Therefore, we obtain the upper bound of A 2 , 
A 2 < l d+1 E Ki < ( 7 d+1 E ^M|n[ 

n n 

l 



(70 
7* 



d+l 



(a-fl) d 



d d xip(a 



(A27) 
oo after 



This tends to when we take the limit I 
7 — !> 0. Hence we can conclude that A 2 —> 0. 

Next, we consider the power-law potential. The power- 
law potential is expressed as K(r) ~ L a ~ d /r a for < 
a < d and its derivative is 



K'{r) 



-d 



r a+l 



We can evaluate KL as 



L a-d L a- 

D1+ l 

n,0 



(al\n 



\a+l ' 



Therefore, when < a < d — 1, 

T ol d / J 



d-a-l 



When a = d — 1 , 



A, < 



— ~ - In ( - 



(A28) 

(A29) 
/ 

I' 

(A30) 
(A31) 



When d — 1 < a < d, since X)n(l/I"1 Q+1 ) is finite in the 
limit L/l —¥ 00, 



Ai<l d+1 J2 K ko^l 



d+ 



a — d 



la+l 



d-c 



(A32) 



In any case, we have A 2 — >• when the conditions (A4) 
are fulfilled. 

In this way, we proved (C7) both for the Kac potential 
and the power-law potential. 



Appendix B: The rigorous justification of the 
saddle-point method 

In Appendix. A, we proved that there exists a function 
g(L, I) such that 



and 



\U-H\< L d g(L,l) 



lim lim g(L, I) = 0. 



(Bl) 
(B2) 



For the definition of H, see Eq. (Al). We evaluate the 
partition function 



2 = EV 



(B3) 



where ^2 denotes the summation over the spin configu- 
ration {<7i} under the restriction of -k at — m. From 
Eq. (Bl), the partition function is bounded by 



E exp 

{m p } 



< z < e ' ex p 

{m„} 



-/3K + E5 , K)-^ d 5(i,0 
p 

-f3H + J2S(m p )+/3L d g(L,l) 



(B4) 



Here, X]{m } denotes the summation over all possible 
values of {m p } with a restriction of ^ m p = m. We 
defined S(m v ) = lnT,, i., cR 1 „ _ m i 1 Let us 



define 



max 
■u 12 P r "p- 



-/3# + E s ( m p) 



Then, a lower bound of Z is given by 
Z > e ~0 F ' -P Ld a{L,i) _ 

An upper bound is also obtained as follows: 

{m p } 
{m p } 



(B5) 



(B6) 



Here, we assume that ^ TO 1 < cl nd with a positive con- 
stant c and a positive integer n, both of which are inde- 
pendent of L and I. We can confirm that this assumption 
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is satisfied in individual models. For example, in the Ising 
model, TO p takes the following value 

l d m p G {l d ,l d -2,l d -4,...,-l d }. 

Therefore, in this case, J2 mp 1 = ^ + 1 < 2/ d . Also in 
other models, we can check this property. Therefore, we 
have 



Z < (d^n e -PF'+l3L d g(L,l) _ 



(B7) 



(note that ft = (L/l) d .) From Eqs. (B5) and (B7), we 
have 

F* , N 1 \n(cl nd ) 11 F* 

(B8) 

Hence, using Eq. (B2), we obtain 

f = —— lim lim — 7 In Z 

P* 

= lim lim — T . (B9) 

Thus, the saddle-point method in Eq. (21) is justified 
rigorously. 



Appendix C: The derivation of the lower bound (23) 

From Eq. (21) and F MF (m,T,H) = H[m] - TS(m), 
we have 



F(m,T,H)-F MF (m,T,H) 

1 H-T [ S(m(x))d d a 
0=m } I J Cd 



mm 

{m(x)|m = 



+ TS{m) 



mm 

{m(x) \rho— m} 



\m r 



+ TS(m). 



1 [ S{m{x))d d i 
Jc d 

(CI) 



Let us define C/ max as the maximum interaction eigen- 
value except for U a = 1, 



U max = maxU n . 
Taking into account of the relations 
^2 U n \m n \ 2 < U max ^ 



n/0 



n#0 



and 



/ m(x) 2 d d x = m 2 + \rh r 



(C2) 



(C3) 



(C4) 



the RHS of (CI) is written as 



\m n \- - T I (S(m(x)) - S(m)) d d x 
>c d 



>^-U max m 2 +TS(m) 



j ■ 

Jc d 



■fc, 



d d x 



-y[/ max m(x) 2 -T5(m(x)) 



F MF (m,T/U max ,H) 
d d xF MF (m(x),T/U max ,H) 



(C5) 



>c d 

Thus, we obtain 

F{m,T,H)>F MF (m,T,H) 

[ d d xF MF (m{x),T/U max ,H) 
Jc d 

(C6) 



+ C^max min 

{m(a:)|mo- m} 



- U mex F MF (m, T/U max , H). 
Here, we can prove 



mm 

{m(x) |mo— m} 



/ d d xF MF {m(x),T/U max ,H) 
Jc d 

= CE{F MF (m, T/(7 max , if)}. (C7) 



Therefore, we obtained the lower bound (23). 

Let us show Eq. (C7). Later, we omit the temperature 
dependence and the external field dependence of the free 
energy. Let us define a function G(m) as 



G(m) 



mm 

{m(x) \rho=m} 



/ d d xF MF (m(x)) 
JC d 



(C8) 



If we consider the uniform conhguration, m(x) = m, we 
have J c d d xFM F (m(x)) = FM F (m). Therefore we ob- 
tain G(m) < FM F (m). From the definition of the convex 
envelope, 



G(m) = min 

{m(x)\mo— m} 

> min 

{m(x)\m — m} 

> min 

{m(x)\mo— m} 



d d xF M F(m(x)) 



c d 
CE 



d d xCE{F M F(m(x))} 
ji^MF (jf d d xm(x) 



= CE{F MF (m)}. 
Hence, we obtain 

CE{F MF (m)} < G(m) < F MF (m). 



(C9) 



(CIO) 



Next let us prove that G(m) is a convex function of m. 
For a constant < A < 1, we consider 



AG(mi) + (1 - A)G(m 2 ). 



(Cll) 
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The functions G{m{) and G(m 2 ) can be written as 

G(mi) = / d d xF MF {m\{x)), / d d xm* 1 {x) — mi, 

(C12) 



G(m 2 ) = / (I xFMFfm^a;)), / d d xm 2 (x) = m 2 , 
Jc d ' Jc d 

(C13) 



where m^a;) and m 2 {x) are the functions that minimize 
the function J c d d xFMF(m(x)) for the fixed global vari- 
ables mi and mi, respectively. Putting \x\ — y\ and 
(1 — A).ti = zi, we obtain 

AG(roi) + (1 - A)G(m 2 ) 
= y rfyi y dx 2 ---dx d F M F (m\ (^-,x 2 , . . . ,x d ^ 

f dZ! [ 
Jo Jo 



+ 



dx 2 ■ ■ ■ dx d F M F ( m* 2 ( - ^ .^2, ■ • • , x d 



= J dxi j dx2---dx d F M ¥(ml(^-,X2,..., 



Xd 



j dx x J dx 2 ■ ■ ■ '/•'•,/ /'mi \^n 2 y - - - 



+ 



x\ — A 



Here, we define a function m*(x) as 
m*(a;) = < 



m* ^^-,a;2, . . . ,x d ) for < x\ < A, 



m 2 ( ~r — ^,£2, . . . ,x d ) for A < xi < 1. 



. 1 - A 

Then, the following condition is satisfied 



Jc d 



xm*(x) — Ami + (1 — A)m-2. 



(C15) 
(C16) 



Therefore, we obtain from (C14) 



AG(mi) + (1 - A)G(m 2 ) = f d d xF MF (m*(x)) 

Jc d 

> G(Xmi + (1 - A)m 2 ). 

(C17) 

We thereby conclude that the function G(m) is convex. 
From Eq. (CIO) and the convexity of G(m), 



G{m) = CE{F MF (m)} 



(C18) 



(C14) holds. 
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